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ABSTRACT 



In the recent work of Nakariakov et al. (2004), it 
has been shown that the time dependences of den- 
sity and velocity in a flaring loop contain pronounced 
quasi-harmonic oscillations associated with the 2nd 
harmonic of a standing slow magnetoacoustic wave. 
That model used a symmetric heating function (heat 
deposition was strictly at the apex). This left out- 
standing questions: A) is the generation of the 2nd 
harmonic a consequence of the fact that the heating 
function was symmetric? B) Would the generation 
of these oscillations occur if we break symmetry? C) 
What is the spectrum of these oscillations? Is it con- 
sistent with a 2nd spatial harmonic? The present 
work (and partly Tsiklauri et al. (20G4b)) attempts 
to answer these important outstanding questions. 
Namely, we investigate the physical nature of these 
oscillations in greater detail: we study their spec- 
trum (using periodogram technique) and how heat 
positioning affects the mode excitation. We found 
that excitation of such oscillations is practically in- 
dependent of location of the heat deposition in the 
loop. Because of the change of the background tem- 
perature and density, the phase shift between the 
density and velocity perturbations is not exactly a 
quarter of the period, it varies along the loop and is 
time dependent, especially in the case of one foot- 
point (asymmetric) heating. We also were able to 
model successfully SUMER oscillations observed in 
hot coronal loops. 

Key words: Sun: flares - Sun: oscillations - Sun: 
Corona - Stars: flare - Stars: oscillations - Stars: 
coronae . 



wind acceleration problems. Observational evidence 
of coronal wa ves and oscillations in EUV are nu- 
merous fe.g. ll(5famn et all Il999t lOfman fc WaneL 
l2n?il lWa,ng et all I2002L l2003D ~ Radio band ob- 
servations also demonstrate various kinds of oscilla- 
tions (e.g., the quas i-periodic pulsations, or QPP, see 
lAschwandenl l|1987t ) for a review), usually with peri- 
ods from a few seconds to tens of seconds. Decimetre 
and microwave observations show much longer peri- 
odiciti es, often in association with a flare. For ex- 
ample, IWang fc Xie (2000) observed QPP with the 
periods of about 50 s at 1.42 and 2 GHz (in as- 
sociation with an M4.4 X-ray flare). Similar pe- 
riodicit ies have been observed in the X-ray band 
(e.g. l|McKenzie fc MiillanL llQQTt iTerekhov et"!!! 
120021) 1 and in the white-hght emission associated 
with the stellar flaring loops ijMathioudakis et al.l 
l2fl03j) . A possible interpretation of these medium 
period QPPs may be in terms of kink or torsional 
modes l|Zaitsev fc Stepanovl ll989l )■ 
In our pre v ious, preliminary study 
([Nakariakov et all |2Q^), we outlined an alter- 
native, simpler, thus more attractive mechanism for 
the generation of long-period QPPs. That model 
used a symmetric heating function (heat deposition 
was strictly at the apex). This left the outstanding 
questions: A) is the generation of the 2nd harmonic 
a consequence of the fact that the heating function 
was symmetric? B) Would the generation of these 
oscillations occur if we break symmetry? C) What 
is the spectrum of these oscillations? Is it consistent 
with a 2nd spatial har mon ic ? The present work 
(and partly Tsiklauri ct^ (|2004bD l attempts to 
answer these important outstanding questions. 

We also were able to model successfully SUMER os- 
cillations observed in hot coronal loops. 



1. INTRODUCTION 



Magnetohydrodynamic (MHD) coronal seismology 
is one of the main reasons for studying waves in 
the solar corona. Such studies also are impor- 
tant in connection with coronal heating and solar 



The paper is organised as follows: in sect. 2 
we present the numerical results, with subsec- 
tion 2.1 dedicated to the case of apex (symmet- 
ric) heating which completes the work started in 
Nakariakov et al. (2004), and subsection 2.2 sum- 
marising our findings in the case of single footpoint 
(asymmetric) heating. In sect. 3 we present some pre- 
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liminary results of modelling of SUMER oscillations 
observed in hot coronal loops. We close with conclu- 
sions in sect. 4. 



2. NUMERICAL RESULTS 



The model that we use to describ e plasma dynam- 
ics in a coronal loop is outlined in ' Nakariakov et al.l 
(2004); Tsiklauri et al. ( 2004a). Here we just add 
that, when numerically solving the ID radiative hy- 
drodynamic equations (infinite magnetic field ap- 
proximation), and using a ID version of the La- 
grangian Re-map code (Arber et al. 2001) with the 
radiative loss li miters, the radiative los s function was 
specified as in iTsikla uri et al"! l l 2004a|) which essen- 
tially is the iR.osner et al.l 1119781) law extended to a 
wider temperature range ijPeres et all Il982t IPriesl 

We have used the sa me heating function as in 
iTsikla^iri et all I^HMi)- The choice of the tempo- 
ral part of the heating function is such that at all 
times there is a small background heating present 
(either at footpoints or the loop apex) which ensures 
that in the absence of flare heating (when a, which 
determines the flare heating amplitude, is zero) the 
average loop temperature stays at 1 MK. For easy 
comparison between the apex and footpoint heating 
cases we fix Qp, flare heating amplitude, at a differ- 
ent value in each case (which ensures that with the 
flare heating on when a = 1 the average loop tem- 
perature peaks at about the observed value of 30 MK 
in both cases). 

In the numerical runs in sect. 2, l/(2cr^) was fixed 
at 0.01 Mm~^, which gives a heat deposition length 
scale, as — 7 Mm. This is a typical value de- 
termined from the observations l)Aschwanden et all 
|200^) . The flare peak time was fixed in all numeri- 
cal simulations at 2200 s. The duration of the flare, 
CTj, was fixed at 333 s. The time step of data visu- 
alisation was chosen to be 0.5 s. The CFL limited 
time-step used in the simulations was 0.034 s. 



2.1. Case of Apex (Symmetric) Heating 

In this subse ction w e complete the analysis started 
inlNakariako v et'aD (j2004), namely for the same nu- 
merical run we study the spectrum of oscillations at 
different spatial points. 

As was pointed out in iNakariakov et alJ l)2004j) . the 
most interesting fact is that we see clear quasi- 
periodic oscillations, especially in the second stage 
(peak of the flare) for the time interval t = 2500 — 
2800 s (cf. Fig. 1 in lNakariakov et all l|2004D 'l. Such 
oscillations are frequently seen during the solar fl ares 
obser ved in X-rays, 8-20 keV, (e.g. Terckhov et^ 
1)2002 1') as we ll as stellar flares ob served in white-light 
(e.g. iMathio^idakis et all l|2003D l. 



VELOCITY PERIODOGRAM (SPECTRUM) 

5x10*1 ' 7> ' I 

4x10+ h '^'(\ H 

3x10* h '/ \\ -J 

2x10*^ 1 H -i 

i /' >\ i 

1x10* h /' \ -4 

0.000 0,010 0,020 0,030 

f [Hz] 

DENSITY PERIODOGRAM (SPECTRUM) 




0,000 0,010 0,020 0.030 

f [Hz] 

Figure 1. Case of apex (symmetric) heating: Peri- 
odogram (spectrum) of the velocity and density oscil- 
latory component times series outputted in the fol- 
lowing three points: loop apex (solid curve), 1/4 
(dash-dotted curve) and 1/6 (dashed curve) of the ef- 
fective loop length (4-8 Mm), i.e. at s = 0,-12,-16 
Mm. 

Before discussing the physical nature of these oscilla- 
tions, it is worth recalling for completeness the sim- 
ple ID analytic theory of standing sound waves. For 
ID, linearised, hydrodynamic equations with con- 
stant unperturbed (zero order) background variables, 
the solutions for density, p, and velocity, Vx, can be 
easily written as 



p{s 




where Cs is the speed of sound, A is wave ampli- 
tude, L is loop length, n = 1,2,3,... is the har- 
monic number, and s is the distance along the loop. 
Note the (relative) phase shift between Vx and p 
is AP/P = -(7r/2)/(27r) = -1/4, where P is the 
standing wave period, while this ratio is zero for 
a propagating w a ve. A lso, Eqs.(3) and (4) from 
INakariakov et all l)2004fl are missing a factor of 2, 
while our Eqs.(l) and (2) correct this previous omis- 
sion. 

In Fig. 1 we present a periodogram (which here we 
use interchangeably with the (power) spectrum, al- 
though strictly speaking the power spectrum is a the- 
oretical quantity defined as an integral over contin- 
uous time, and of which the periodogram is simply 
an estimate based on a finite amount of discrete data 
(cf. IScargld l)1982 t) and his Eq.(lO) in particular)) of 
the velocity and density time series outputted at the 
three points: loop apex, 1/4 and 1/6 of the effective 
loop length (48 Mm), i.e. at s = 0,-12,-16 Mm. 
The first two points are chosen to test whether the 
simple analytic solution for ID standing sound waves 
(see below) is relevant in this case. The third point 
(1/6) was chosen arbitrarily (any spatial point along 
the loop where density and velocity of the standing 
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Figure 2. Case of apex (symmetric) heating: os- 
cillatory components of time series outputted at ±6 
Mm and ±18 Mm in the time interval of 2500-2800 
s. The solid curve shows plasma number density in 
units o/lO^^ cm~^ . The dashed curve shows velocity 
normalised to 400 km s~^ . 

waves does not have a node would be equally ac- 
ceptable) . As expected for a 2nd spatial harmonic of 
a standing sound wave in the velocity periodogram 
there are two clearly defined peaks and the largest 
peak corresponds to 1/4 of the effective loop length, 
while the smaller peak corresponds to 1/6. Note that 
at the loop apex the periodogram gives (solid line is 
too close to zero to be seen in the plot). The density 
periodogram shows the opposite behaviour to that of 
the velocity with the largest peak corresponding to 
the loop apex, while 1/6 of the effective loop length 
corresponds to a smaller peak, and 1/4 is close to 
zero. The locations of the peaks are at about 0.0155 
Hz i.e. the period of the oscillation is 64 s. The pe- 
riod of a 2nd spatial harmonic of a standing sound 
wave should be 

P = L/C, = i/(l. 52 X 10^\/r), (3) 

where T is plasma temperature measured in MK, 
while L is in meters. If we substitu te an effective loop 
lengt h L = 48 Mm (see Fig. 2 in ijNakariakov et all 
l200l 'l and an average temper ature of 25 M K (see 
top panel in Fig. 1 inlNakariak ov et "all (|2004l) in the 
range of 2500-2800 s - the quasi periodic oscillations 
time interval we study) we obtain 63 s, which is close 
to the result of our numerical simulation. Such a 
close coincidence is surprising bearing in mind that 
the theory does not take into account variation of 
background densi ty and velocity oy e r time , while we 
see from Fig. 1 in lNakariakov et'aLl l)2004l) that even 
within a short interval of a flare, i.e. 2500-2800 s, all 
physical quantities vary significantly with time. To 
close our investigation of the physical nature of the 
oscillations we study the phase shift between the ve- 
locity and density oscillations and compare our simu- 
lation results with analytic theory. In Fig. 2 we plot 
time series, outputted at ±6 and ±18 Mm, of the 
plasma number density in units of 10^^ cm~^ and 
velocity, normalised to 400 km s~^. These points 
were selected so that one symmetric (with respect to 



the apex) pair (±6 Mm) is close to the apex, while 
another pair (±18 Mm) is closer to the footpoints. 
We choose these pairs because we wanted to compare 
how the phase shift is affected by spatial location. 
One would expect stronger upflows close to foot- 
points (due to chromospheric evaporation) , which in 
turn alters the phase shift. Note that phase shift be- 
tween the density and velocity is different (see below) 
for standing and propagating (with flows) acoustic 
waves). We gather from the graph that: (A) clear 
quasi-periodic oscillations are present; (B) they are 
shifted with respect to each other in time; (C) near 
the apex (±6 Mm) the phase shift is close to that 
predicted by ID analytic theory; (D) close to the 
footpoints (±18 Mm) the phase shift is somewhat dif- 
ferent from the one predicted by ID analytic theory. 
In the last case the discrepancy can be attributed to 
the presence of flows near the footpoints. The main 
reason for the overall deviation is due to the fact that 
analytic theory does not take into account variations 
of background density and velocity in time and that 
density gradients in the transition region are not pro- 
viding perfect reflecting boundary conditions for the 
formation of standing sound waves. 

Another interesting result is that even with the wide 
variation of the parameter space of the flare, its dura- 
tion, peak average temperature, etc., we always ob- 
tained a dominant 2nd spatial harmonic of a stand- 
ing sound wave with some small admixture of 4th 
and sometimes 6th harmonics. Our initial guess was 
that this is due to the symmetric excitation of these 
oscillations (recall that we use apex heat deposition). 
In order to investigate the issue of excitation further 
we decided to break the symmetry and put the heat- 
ing source at one footpoint, hoping to see excitation 
of odd harmonics 1st, 3rd, etc. 



2.2. Case of Single Footpoint (Asymmetric) 
Heating 

For single footpoint heating we flx sq = 30 Mm in 
Eq.(l) in Nakari akov et al. (2Q^), i.e. (spatial) peak 
of the heating is chosen to be at the bottom of the 
transition region (top of chromosphere). Initially we 
run a code without flare heating, i.e. we put a = 
(in this manner we turn off flare heating). Eq = 0.02 
erg cm~'^ was chosen such that in the steady 
(non-flaring) case the average loop temperature stays 
at about 1 MK. Then, we run the code with flare 
heating, i.e. we put a = 1, and flx Qp at 1 x 10^, so 
that it yields a peak average temperature of about 
30 MK. The results are presented in Fig. 3. Dur- 
ing the flare the apex temperature peaks at 38.38 
MK while the number density at the apex peaks at 
3.11 X 10"'^^ cm~^. In this case, as opposed to the case 
of symmetric (apex) heating, the velocity dynamics 
is quite different. Since the symmetry of heating is 
broken there is a non-zero net flow through the apex 
at all times. However, as in the symmetric heating 
case, we again see quasi-periodic oscillations super- 
imposed on the dynamics of all physical quantities 
(cf time interval of t = 2400 - 2700 in Fig. 3). 
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Figure 3. Case of single footpoint ( asymmetric ) heat- 
ing: Average temperature, temperature at apex, and 
number density at the apex as functions of time. 

In Fig. 4 we present time-distance plots of veloc- 
ity and density for the time interval 2400-2700s, 
where the quasi-periodic oscillations are most clearly 
seen. Here we again subtracted the slowly vary- 
ing background (with respect to oscillation period). 
The picture is quite different from the case of apex 

S symmetric) heating (c ompare it with Fig. 2 in 
Nakaria kov et al.ll2004() \ This is because now the 
node in the velocity (at the apex) moves back and 
forth periodically along the apex, and at later times 
{t > 2550 s) stronger flows are now present. How- 
ever, the physical nature of the oscillations remains 
mainly the same. i.e. a 2nd spatial harmonic of a 
sound wave, but now with an oscillating node at the 
apex. 

To investigate this further we plot in Fig. 5 a peri- 
odogram (spectrum) of the velocity and density os- 
cillatory component time series outputted at the fol- 
lowing three points: loop apex, 1/6 and 1/4 of the 
effective loop length. We gather from the graph that 
the periodogram (spectrum) is more complex than in 
the case of apex (symmetric) heating. In the veloc- 
ity periodogram at the apex there is a peak with a 
frequency higher than that of 2nd spatial harmonic 
of a standing sound wave. This is the frequency with 
which the node of the velocity oscillates (see discus- 
sion in the previous paragraph). It has nothing to 
do with the standing mode, but is dictated by the 
excitation conditions of the loop which acts as a dy- 
namic resonator. Let us analyse now how this pe- 
riodogram compares with ID analytic theory. The 
peak in the periodogram corresponding to 1 /6 of the 
effective loop length (dashed line) corresponds to a 
frequency of about 0.017 Hz, i.e. the period of os- 
cillation is 59 s. Again, the period of a 2nd spatial 
harmonic of a standing sound wave should be cal- 
culated from Eq. (3). If we substitute the effective 
loop length L = 48 Mm (see Fig. 4) and an average 
temperature of 26 MK (see top panel in Fig. 3 in the 
range of 2400-2700 s) we obtain 62 s, which is close 
to the result of our numerical simulation (59 s). 

Next, we studied the phase shift between velocity 



Density 




Distaoce [Mm] Distance [Mm] 

Figure 4- Time-distance plots of the velocity and 
density oscillatory components in the time interval 
of 2400-2700s for the case of single footpoint (asym- 
metric) heating. 
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Figure 5. As is Fig. 1 but for the case of single foot- 
point (asymmetric) heating. Time interval here is 
2400- 2700s. 



and density oscillations, and compare our simulation 
results with the ID analytic theory. In Fig. 6 we pro- 
duce a plot similar to Fig. 2, but for the case of asym- 
metric heating. The deviation, which is greater than 
in the case of apex (symmetric) heating, can again 
be attributed to the over-simplification of the ID an- 
alytic theory, which does not take into account time 
variation of the background physical quantities and 
imperfection of the reflecting boundary conditions 
(see above). More importantly, in the asymmetric 
case strong flows are present throughout the flare 
simulation time. Thus, if linear time dependence is 
assumed, which is relevant within the short interval 
2400-2700 s of the flare, then Eqs.(l)-(2) would be 
modified such that phase shifts would vary secularly 
in time. This is similar to that seen in Fig. 6. 

Yet another interesting observation comes from the 
following argument: in a steady ID case analytic the- 
ory predicts that the phase shift between the density 
and velocity should be (A) zero for for a propagat- 
ing acoustic wave and (B) quarter of a period for 
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Figure 6. As in Fig. 2, but for the case of single 
footpoint (asymmetric) heating. Time interval here 
is 2400-2700S. 



Figure 8. As in Fig. 7, but with with subtracted 
slowly varying background. 
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Figure 7. Case of apex (symmetric) heating: Aver- 
age temperature, temperature at apex, number den- 
sity at the apex, and velocity at the apex as functions 
of time. 



the standing acoustic wave. Since in the asymmet- 
ric case strong flows are present, we see less phase 
shift between the velocity and density in Fig. 6 as 
one would expect. 

Thus, the resu l ts of th e present study (and partly 
iTsiklauri et al.l l)2004bl )'l provide further, and more 
definitive proof than in Nakariakov ct al. ( 20Q4) that 
these oscillations are indeed the 2nd spatial harmonic 
of a standing sound wave. However, the present work 
also reveals that in the case of single footpoint (asym- 
metric) heating the physical nature of the oscillations 
is more complex, as the node in the velocity oscillates 
along the apex and net flows are present. 



SUMER OSCILLATIONS OF HOT 
CORONAL LOOPS 



Recently hot (T = 6.3 MK and T = 8 MK) coro- 
nal loop osci llations were obse rved with SUMER (cf. 
for review iWang et all l|2n0.lM These were inter- 



preted as standing acoustic waves lOfman fc Wand 
l)2002f) . So far only highly simplified models (isother- 
mal, no stratification, no Helium presence, no ex- 
ternal heating input, no transition region or chro- 
mosphere, etc.) were used to study these oscilla- 
tions. Excitation of the standing waves was done 
artificially perturbing velocity of plasma, globally, 
along the loop and then their evolution (rapid de- 
cay) was modelled numerically with heat conduction 
being dominant factor Ofman & Wang ( 2002) . Here 
we use much more sophisticated model, as described 
above, to self-consistently excite these oscillations by 
adding external heat input. In a similar manner as 
in previous sections we use apex (symmetric) heat- 
ing. Parameters of the heating function were chosen 
such the average temperature peaks at about 8 MK. 
The time dynamics of various quantities are shown 
in Fig. 7. In Fig. 8 we plot the same quantities 
with subtracted slowly varying background. This 
shows clear, periodic, decaying oscillations are ex- 
cited. Loop length was fixed at L — 350 Mm. 

From Fig. 8 it follows that the period of oscillations 
is about 800 s « 13 min, which is similar to the 
observed range of 17.6 ± 5.4 min SU MER oscilla- 
tions (based on 54 Doppler-shift cases ijWang et al.l 
I2003M . Moreover, the same result is obtained from 
the simple ID analytic theory when loop length 
L ~ 350 Mm and temperature of 8 MK is substi- 
tuted into Eq.(3). Thus, our model seems to describe 
adequately essential observed features of these oscil- 
lations. Parametric study showed that in the case 
SUMER oscillations heating function should have 
fairly short temporal width (about 1 minute). The 
main difference from previous models is that we self- 
consistently deal with the excitation problem. 

Note also that in this simulation period of the oscilla- 
tions is increased. This can be attributed to two fac- 
tors: (A) decrease of backround tempreature (which 
would not be observed by SUMER because it looks 
at lines with fixed temperatures (T = 6.3 MK and 
T — 8 MK)), and (B) dissipative effects (mainly heat 
conduction) . It seems that the increase in period ob- 
served by SUMER should be attributed to the latter 
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effect (B). 



4. CONCLUSIONS 

Initially we used a ID radiative hydrodynamics 
loop model which incorporates the effects of grav- 
itational stratification, heat conduction, radiative 
losses, added external heat input, presence of he- 
lium, hydrodynamic non-linearity, and bulk Bra- 
ginskii viscosity to simulate flares ijTsiklauri et all 
I^HmH)- As a byproduct of that study, in prac- 
tically all the numerical runs quasi-periodic os- 
cillations in all physica l quantities were detected 
ijNakariakov et all l2004j) . Such oscillations are fre- 
quently seen during th e solar flares obser ved in 
X-rays, 8-20 keV (e.g. iTerekhov etlll ^2)) as 
well as stell ar flare s obse rved in white-light (e.g. 
iMathioudakis et al.l i200.lD1. Our pr esent analy- 
sis (and partlv iTsiklauri et alJ l)2004b() ') shows that 
quasi-periodic oscillations seen in our numerical sim- 
ulations bear many similar features compared to ob- 
servational datasets. In this work we tried to an- 
swer important outstanding questions (cf. Introduc- 
tion section) that arose from the previous analysis 
l|Nakariakov et all l200^ . 

In sumr nary the present study fand'Nakariak ov et al.l 
(200l): ITsiklauri et alJ p004b) ) established the fol- 
lowing features: 

• We show that the time dependences of density 
and temperature in a flaring loops contain well- 
pronounced quasi-harmonic oscillations associ- 
ated with standing slow magnetoacoustic modes 
of the loop. 

• For a wide range of physical parameters, the 
dominant mode is the second spatial harmonic, 
with a velocity oscillation node and the density 
oscillation maximum at the loop apex. This re- 
sult is practically independent of the positioning 
of the heat deposition in the loop. 

• Because of the change of the background tem- 
perature and density, and the fact that den- 
sity gradients in the transition region are not 
providing perfect reflecting boundary conditions 
for the formation of standing sound waves, the 
phase shift between the density and velocity per- 
turbations is not exactly equal to a quarter of a 
period. 

• We conclude that the oscillations in the white 
light, radio and X-ray light curves observed dur- 
ing solar and stellar flares may be produced 
by slow standing modes, with the period deter- 
mined by the loop temperature and length. 

• For a typical solar flaring loop the period of os- 
cillations is shown to be about a few minutes, 
while amplitudes are typically of a few percent. 

• Our model seems to describe adequately essen- 
tial observed features of SUMER oscillations. 



Parametric study showed that in this case heat- 
ing function should have fairly short temporal 
width (about 1 minute). 

The novelty of this study is that by studying the 
spectrum and phase shift of these oscillations we pro- 
vide more definite proof that these oscillations are in- 
deed the 2nd harmonic of a standing sound wave, and 
that the single footpoint (asymmetric) heat position- 
ing still produces 2nd spatial harmonics, although it 
is more complex than the apex (symmetric) heating 
as the node in the velocity oscillates along the apex 
and net flows are also present. 
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